Singular non-ordering susceptibility at a Pomeranchuk instability 
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We study magnetic susceptibilities of two-dimensional itinerant electron systems exhibiting symmetry- 
breaking Fermi surface distortions, the so-called rf-wave Pomeranchuk instability, in a magnetic field. In a 
pure forward scattering model, the longitudinal susceptibility x~ is found to exhibit a jump at a critical point. 
The magnitude of this jump diverges at a tricritical point. When scattering processes involving finite momen- 
tum transfers are allowed for, x~ is expected to diverge also at a critical point. The system displays multiple 
critical fluctuations. We argue that the features of x~ are general properties associated with singularities of a 
non-ordering susceptibility, leading to implications for a variety of materials including Sr 3 Ru 2 7 . 

PACS numbers: 75.40.-s, 71.18.+V, 05.50.+q, 64.60.F- 



I. INTRODUCTION 

The electronic nematic is currently among the most intrigu- 
ing and intensely investigated states of quantum matter. Un- 
like its classical counterpart [1] this phase originates from 
interactions between structureless, point-like particles (elec- 
trons). On the other hand, just like in the case of classical 
nematics, these systems exhibit orientational symmetry break- 
ing while preserving translational invariance. The electronic 
nematic state is reached via a phase transition from either the 
electronic smectic state corresponding to the so-called charge- 
stripe order [2 or the uniform state as found in the t-J |3[] and 
Hubbard Si models. The latter route is accompanied by 
spontaneous Fermi surface symmetry breaking with a c/-wave 
order parameter, and is often referred to as the cZ-wave Pomer- 
anchuk instability Ja] . 

The occurrence of the electronic nematic was indicated in 
a number of correlated electron systems. In two-dimensional 
electron gases, the electronic nematic is signaled by the strong 
anisotropy of magnetotransport at low temperature in half- 
filled higher Landau levels 0]. For cuprate superconductors, 
there exists a strong tendency toward the nematic state despite 
competition with cooper pairing formation, as seen in the pro- 
nounced anisotropy of magnetic excitation spectra @]. The 
bilayer ruthenate Sr3Ru207 displays the evidence of the ne- 
matic state II9[ llOll . which is realized in a certain range of a 
magnetic field bounded by first order transition lines at low 
temperature and a second order transition line at high temper- 
ature; the end points of the second order transition line are 
tricritical points ll ill . 

Motivated by experiments performed for Si"3Ru207, we 
consider the nematic instability from a uniform state, namely 
the Pomeranchuk instability (PI), in the presence of a mag- 
netic field. We find that the PI is accompanied by peculiar 
singularities of the uniform spin susceptibility, which can be 
observed in real materials like S^Ri^O?. 

Our results are not specific to the PI and can be understood 
as universal singular behavior of a non-ordering susceptibil- 
ity. The concept of a non-ordering susceptibility appears in 
the theory of tricriticality developed for helium mixtures and 
some metamagnetic materials fl2l[l3ll . At the tricritical point 
three different phases become identical lfl4ll and the full phase 



diagram is described by temperature (T), a control parameter 
called a non-ordering field, and an ordering field which is con- 
jugate to the order parameter [15]. The framework naturally 
involves two different kinds of susceptibilities, the ordering 
and non-ordering ones. The former is given by the second 
derivative of the free energy with respect to the ordering field, 
whereas the latter is the second derivative with respect to the 
non-ordering field; the non-ordering susceptibility is thus de- 
fined uniquely in the system in question. In the context of 
Sr^RuaOv, the ordering susceptibility corresponds to the d- 
wave charge compressibility II 1 611 and the non-ordering one to 
the longitudinal magnetic susceptibility. 

In this paper, we first study the uniform magnetic suscep- 
tibilities in a mean-field model of the PI on a square lattice 
in the presence of a magnetic field. We show that the lon- 
gitudinal susceptibility x zz exhibits a jump at a critical point 
(CP) and a power-law divergence at a tricritical point (TCP), 
while the transverse susceptibility x ± exhibits a cusp. In two- 
dimensional systems, the mean-field theory breaks down in a 
critical regime. We then argue that beyond the mean-field the- 
ory x zz diverges at both CP and TCP, which is associated with 
multiple critical fluctuations, i.e., soft Fermi surface fluctua- 
tions and magnetic fluctuations. The singular behavior of x zz 
is interpreted as a general feature associated with singularities 
of a non-ordering susceptibility at a phase transition, leading 
to implications for a variety of materials including Sr^RmOv. 



II. MODEL AND FORMALISM 



We consider the following Hamiltonian for the PI lfl7i [l8ll 
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describing tight-binding electrons on a square lattice, inter- 
acting via a forward scattering interaction and subject to an 
external magnetic field h. The dispersion e£ = -2[?(cos£ x + 
cos k y ) + 2?' cos k x cos k y ] retains the nearest and next-nearest 
neighbor hopping amplitudes; yt is the chemical potential; 
is the electron number operator with momentum k and spin <x; 



N is the total number of lattice sites; the effective interaction 
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involves the li-wave form factor d k - cos k x -cos k y and drives 
the PI. This interaction mimics the actual structure of the ef- 
fective charge interactions in the forward scattering channel 
as obtained within more fundamental models JH-Hl]. The or- 
der parameter corresponding to the PI is given by 
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X zz = 
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FIG. 2: Graphical representation of X" and x ± f° r the Hamiltonian 
(T). The vertex with a circle (square) indicates the form factor 1 (rf k ). 
The solid line denotes electron propagator with the dispersion <fko-- 



The transition can be tuned along an isotherm either by vary- 
ing the chemical potential fi or the magnetic field /; in the 
Hamiltonian (Q}. The latter mimics the actual experimental 
setup in SrjRuaO?. In fact, the model (1) and different ver- 
sions of it turned out to capture numerous aspects possibly 
associated with the Si"3Ru207 compound II 17142611 . We fix the 
band parameters t' = 0.35, g = 1, and y, = 1 in units set by 
t lfl7lTT8ll . The theoretical phase diagram is shown in Fig. [T] 
The PI is realized in a region around the van Hove energy 
(h = 0.8) of the up-spin band. A second-order phase transi- 
tion occurs at high T, but terminates at a TCP, below which 
the transition becomes first order. 
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FIG. 1: The phase diagram in the plane of h and T for f - 0.35, 
g = 1, and ix = 1. T c (7^ st ) denotes a second (first) order transition 
and T a i is a tricritical point. The dashed line represents the position 
of the van Hove energy. 



We now analyze the uniform longitudinal magnetic suscep- 
tibility^ 4 " within the model ([TJ. The standard linear response 
theory yields 
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Here <■ • • ) denotes the equilibrium expectation value with the 
Hamiltonian ([1), [■, ■] is the commutator; 6 a positive infinites- 
imal; S(q, t) = e lH 'S(q)e~ lHt and S(q) is the spin operator in 
momentum space. In the random phase approximation (RPA), 
X zz is given graphically by Fig.|2]and calculated as 
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where £ kcr = e£ — fi + rjd k - ? h; f(x) is the Fermi-Dirac distri- 
bution and its first derivative is denoted by /'. The presence 
of the second term in Eq. (0 is crucial to the present study 
and was overlooked in Ref. [l7l On the other hand, the uni- 
form transverse susceptibility is expressed by a simple bubble 
(Fig.|2) in the RPA, that is, 
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dfe-*<[S + (q, t),S~(-q, 0)]) (6) 



N 

2m 
~~h 



(7) 
(8) 



where S ± = S x ± iS y and m = ^ Ek,<r vfUhur) is the magne- 
tization. 



m. RESULTS 

We plot the susceptibility^ versus temperature in Fig. [3 a) 
as the second-order transition line is crossed at fixed h. The 
quantity x zz exhibits a discontinuity at T = T c (h). The mag- 
nitude of the jump diverges as the field h is varied toward the 
tricritical value hoi = 0.983 [Fig. Ob)]. It follows the power 
law^ zz oc \T — T^r 12711 . The exponent 1/2 is visible only 
sufficiently close to the TCP and very high numerical accu- 
racy is required. The divergence of x zz is observed only in the 
symmetry -broken phase, while in the disordered phase it re- 
mains finite and smooth. In Fig. 0c), the jump of x zz (Ax zr ) 
is plotted along the transition line (solid line). The magni- 
tude of the jump is strongly suppressed around h where T c 
becomes maximal (see also Fig.[T) and is gradually enhanced 
away from it. While x zz necessarily shows a jump upon cross- 
ing a first order transition line, we also plot its magnitude in 
Fig. [3fc) with open circles. It is interesting to note that the 
magnitude of the jump of x zz at a CP is typically comparable 
to that at a first order transition. In contrast to the results for 
X zz , the transverse susceptibility x ± [Fig. HJd)] exhibits a cusp 
upon crossing the T c line at both CP and TCP. Although x ± 
can be enhanced below T c for a different choice of parame- 
ters, its non-analyticity at T c is always a cusp. 

As shown in Fig. [2] x zz involves contributions from two 
diagrams. Its peculiar behavior originates from the second di- 
agram, which is evaluated as the second term in Eq. ([5]). This 
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FIG. 3: (a) and (b) The longitudinal susceptibility^" plotted versus temperature as the second order transition is approached at h = 0.88 (a) 
and at the tricritical field h tIl = 0.983 (b). The dashed line is the result for g=0. (c) The jump ofx lz (A^ z ) along a second order transition line 
(solid line) and at a first order transition (circle) on the side of the symmetry-broken phase. The magnitude of the jump is scaled by x zz on the 
side of the normal phase (x z , z )- (d)^* is plotted for h = 0.88 and 0.983 in a similar fashion to (a). 



term can become finite only in the symmetry-broken phase 
because of the af-wave form factor in the numerator. Hence 
the singular behavior of x zz 1S found only in the symmetry- 
broken phase [Fig.(3ja) and (b)]. This feature was observed 
also in different systems involving non-ordering susceptibili- 
ties considered in the context of helium mixtures l28ll . On the 
other hand, the forward scattering interaction does not con- 
tribute to^ (Fig. [2j, yielding behavior very distinct from^-' 
[Fig. EM)]. 

A. Comparison with a jump of the specific heat 

It is instructive to compare the jump of x zz with a jump 
of the specific heat coefficient (C/T) 11711 . which is plotted 
along the transition line with a solid line in Fig. |4] The jump 
is pronounced around the van Hove energy (h = 0.8) and is 
suppressed monotonously away from it. The quantity C/T 
rapidly increases and ultimately diverges upon approaching 
the vicinity of the TCP; the singularity is given by C/T ~ 
\T - r tri r 1/2 at h = Ki or C/T ~ \h - h^\~ 1!2 at T = T ai 
sufficiently close to the TCP, displaying the same exponent as 
in the case of^ :z . A crucial difference from Fig.[3jc) is that 
while the specific heat exhibits the same singularities even for 
ft = if the PI is tuned, e.g., by the chemical potential iflrjll . 
the jump of x zz at a CP and its divergence at a TCP occur 
only for h + 0. This feature can be checked by observing 
that the second term in Eq. (f5]) vanishes for h = 0. Therefore 
the singular behavior of x zz is a field-induced anomaly. As we 



shall show below, it is a manifestation of a general feature of a 
non-ordering susceptibility and hence the presence of a finite 
non-ordering field (a magnetic field here) is necessary. 
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FIG. 4: The jump of C (AC) along a second order transition line 
(solid line) and at a first order transition (circle) on the side of the 
symmetry-broken phase. The magnitude of the jump is scaled by C 
on the side of the normal phase (C„). 



B. Scaling analysis and criticality beyond the RPA 

So far we have considered the forward scattering model. Its 
universal features in the vicinity of CP and TCP are well de- 
scribed by mean-field theory, and are correctly captured also 
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by a Landau-type </> 6 theory |[L2tl . By allowing a small mo- 
mentum transfer in Eq. ©, we can incorporate order parame- 
ter fluctuations. The obtained singularities of x zz and C/T are 
expected to be modified in a critical regime. By virtue of uni- 
versality, the character of the non-analyticities and the exact 
values of critical exponents can be deduced without perform- 
ing explicit calculations. The order parameter describing the 
PI [Eq. ©] is characterized by the Ising symmetry. The ex- 
act critical exponents of the Ising universality class are well 
known for ordering quantities. The clue to obtain the exact 
exponents of non-ordering quantities is seen in theory of tri- 
critical points 01211 . The scaling hypothesis allows one to write 
the singular part of the Gibbs potential around a TCP as 



G smg (t,g,s) = \g\ 2 - a g(t/\gf,s/\ g \ A ), 



(9) 



where t — (T - T^/T^ and s is an ordering field which is 
conjugate to the order parameter; the scaling field g is chosen 
as g — (h—htn)-at where h is a non-ordering field and <T X is a 
slope of the transition line at the TCP. We focus on a situation 
where the transition line is approached with a finite angle. 

Singular contributions to non-ordering quantities are ob- 
tained by taking derivatives of G S me with respect to hath = h trl 
and 5 = with t — > 0: 



with /?2 = 1 - o and 



X 



dG sin g 

dh 



d 2 G siD 
dh 2 



\tf- 



(10) 



(11) 



with y2 = oc, where the additional subscript "2" indicates the 
exponents corresponding to non-ordering quantities. In the 
present case, the non-ordering field h corresponds to a mag- 
netic field and thus m (x) is the magnetization (longitudinal 
magnetic susceptibility). Since the specific heat coefficient 
C/T at h = /%; is given by 



C/T = — 



9 Gsin; 

dT 2 



(12) 



j2 is the same as the specific heat exponent. Hence we obtain 
the following scaling laws: 



a + 2#> + 72 = 2, 



72 ■ 



(13) 



For an ordinary CP, replacing a TCP with a CP and employing 
the standard form 



G»n g (g,s) = \g\ 2 - a g{s/\gf), 



(14) 



we obtain the same scaling laws. Therefore Eq. ([TBI holds 
for both CP and TCP. This is worth emphasizing that even the 
existence of TCP is not necessary for the fulfillment of the 
relation Eq. ([T3l . 

Since the exact value of a is known both at CP and TCP 
1H2I1 . Eq. ([T"3l gives the exact values of fa an d 72 (Table |I). 
That is, x zz exhibits a divergence in a critical regime at both 
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TABLE I: Critical exponents at a CP and a TCP. Mean-field values at 
a TCP are exponents on the symmetry-broken side. The exponents of 
ordering quantities are also listed for completeness: rj ~ \tf and K d ~ 
|f|~ y , where in the present case r] corresponds to the order parameter 
of the PI [Eq. ([3}] and K d to the <i-wave compressibility [16]. 



CP and TCP, indicating that the PI in a magnetic field is ac- 
companied not only by soft Fermi surface fluctuations but also 
by sizable ferromagnetic fluctuations. It should be noted that 
the singularity of ^ seen in Fig. [3(d) is expected to remain 
the same even in the presence of fluctuations, because the re- 
lation x ± — 2m /h [Eq. ([H)] still holds beyond the RPA as long 
as the effective interaction between electrons retains the spin- 
rotational invariance [29], and m/h is continuous at both CP 
and TCP. 



IV. DISCUSSION 

A. Relevance to Sr 3 Ru 2 07 

There is growing evidence that the bilayer ruthenate 
Sr 3 Ru 2 7 exhibits the PI llgj-TTTh. Although other interactions, 
not included in our model (1), may exist in Sr3Ru207, by 
virtue of universality, the present approach is sufficient to cor- 
rectly capture the universal singularities (Table associated 
with the PI. Inclusion of other terms like spin-orbit coupling 
does not influence the order parameter symmetry [Eq. ([3j] 
and therefore the system can be presumed to remain in the 
two-dimensional Ising universality class even if the micro- 
scopic model is tailored to capture more specific features of 
Sr 3 Ru 2 7 . 

The specific heat jump was reported recently ifloll at a CP 
of the PI. The actual critical region seems tiny enough not to 
be resolved in the experiment and the observed jump corre- 
sponds to the mean-field behavior. Therefore we predict the 
same behavior also for^ [Figs. 3 (a) and (c)]. To the best 
of our knowledge, St"3Ru207 can provide the first example of 
a jump of the magnetic susceptibility at a continuous phase 
transition in condensed matter. On the other hand, at a TCP 
we predict a divergence of the jump ofx zz [Figs. 3(b) and (c)] 
as well as that of C/T (Fig. |4)l, although the latter becomes 
visible only sufficiently close to a TCP, which requires high 
accuracy of an experimental measurement. The quantity x ± 
is not a non-ordering susceptibility and its singularity might 
depend on microscopic details of the system. In particular, 
in Sr3Ru2C<7 the spin-rotational invariance is broken due to 
the spin-orbit coupling. Nonetheless, as far as the underlying 
interaction driving the PI [Eq. (0] does not contribute to x ± 
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(Fig. |2]), we expect a pronounced difference between^ and 
X ± at the PI. 

B. Implications for other materials 

As seen from the derivation of Eq. (jT3j, the singular behav- 
ior of x zz is not a specific feature of the PI in the presence of a 
magnetic field, but can be interpreted as a phenomenon asso- 
ciated with singularities of a non-ordering susceptibility at a 
phase transition in general. This observation provides impli- 
cations for a variety of materials. In particular, when a critical 
region is accessed by experiments, the mean-field predictions 
break down. In such a region, both ordering and non-ordering 
susceptibilities diverge simultaneously (y and j2 in Table H). 
The system features multiple critical fluctuations. The effect 
around a TCP may be more pronounced because the value of 
ji is close to y. Though the universality class is different from 
the present two-dimensional Ising type, singular behavior of a 
non-ordering susceptibility was indeed measured at a TCP in 
some metamagnetic compounds such as CsCoCl3-2D20 and 
DyaAlsOn in 1970s [13]. A manifestation of multiple singu- 
lar susceptibilities is seen also in recent literature. For exam- 
ple in Ref. |30j peculiar behavior of the uniform magnetic sus- 
ceptibility was invoked to propose an explanation of anoma- 
lous properties of YbRh2Si2 and other antiferromagnetic com- 
pounds. Besides these materials, singular behavior of a non- 
ordering susceptibility, which has not been considered much 
in the modern condensed matter field, may provide an inter- 
esting route to interpret anomalous behavior observed in cor- 



related electron systems. 



V. SUMMARY 

We have studied the uniform magnetic susceptibility at the 
PI on a square lattice in the presence of a magnetic field. In 
a mean-field model, x zz is found to exhibit a jump at a crit- 
ical point and a power-law divergence at a tricritical point, 
whereas y ± exhibits a cusp. The mean-field results may be 
modified in a critical region. In particular, x zz is expected 
to diverge also at a critical point, indicating that the PI in a 
magnetic field is accompanied not only by soft Fermi surface 
fluctuations but also by sizable ferromagnetic fluctuations. In- 
voking the scaling hypothesis we have argued that the features 
of are not specific to the PI in a magnetic field, but are gen- 
eral properties associated with singularities of a non-ordering 
susceptibility at a phase transition. While such properties have 
not been recognized well in modern condensed matter litera- 
ture, it can be tested for S^R^Ch by measuring the longitu- 
dinal magnetic susceptibility. 
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